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O ! Abstract 

We consider the classical equations of the Born-Infeld-Abelian-Higgs 
model (with and without coupling to gravity) in an axially symmetric 
D . ansatz. A numerical analysis of the equations reveals that the (gravi- 



tating) Nielsen-Olesen vortices are smoothly deformed by the Born-Infeld 



^ . interaction, characterized by a coupling constant (3 , and that these so- 

Cd I lutions cease to exist at a critical value of (3 . When the critical value is 

approached, the length of the magnetic field on the symmetry axis becomes 
infinite. 



1 Introduction 



Recently, various gauge field theories were studied in which the dynamics of the 
gauge fields is governed by a Born-Infeld (BI) lagrangian The relationship 
between string theory and BI lagrangian obtained in started an important 
activity in the topic which is reviewed e.g. in |P, §]. 

Among the various aspects of BI field theory, the construction of classical 
solutions occupies a large place. Many investigations were carried out in Abelian 
as well as in non- Abelian [^, ^ models. Often, the pattern of classical solutions 
is enriched by the occurence of the Bl-term. The case of the SU(2) gauge theory 
is typical of this type: no static finite energy solutions exist with the minimal 
Yang-Mills action but solutions exist when the BI term is included 0. Lately, 
various Born-Infeld solitons have been considered in the context of gravity (see. 
e.g. Bi). 

Perhaps the simplest model containing solitons for which the effects of the BI 



term may be studied is the Abelian-Higgs model [|10| in which vortex-solutions 
are available. Several aspects of Bl-vortex solutions have been considered so 



far, namely the construction of self dual equations from the BI lagrangian [|TT 



and the construction of electrically charged vortex solutions . Very interesting 
applications of the vortex solutions are cosmic strings ||T^ which are believed to 
have played a role in the early phase of the Universe; it is therefore natural to 
also consider the effect of gravitation on the Bl-vortex solutions. 

Here we consider the Abelian-Higgs-Born-Infeld model coupled to gravity and 
we preform a detailed numerical analysis of the equations in the axially symmet- 
ric ansatz. We find that, when the various coupling constants of the Einstein- 
Abelian-Higgs equations are such that a magnetically charged soliton exists (i.e. 
for (5"^ = oo where (3"^ denotes the BI coupling constant), then the solution is 
smoothly deformed by the BI interaction. When the parameter approaches 
a critical value, the length of the magnetic field of the solution on its symmetry 
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axis becomes infinite and the solution disappears. 

The paper is organized as follows : in Sect. 2 we describe the Lagrangian, 
the ansatz and the classical equations. Several known solutions and some of their 
features are summarized in Sect. 3. Sect. 4 contains our results on the evolution 
of these solutions in presence of the Born-Infeld interaction. 



2 The equations. 

We consider the gravitating Born-Infeld- Abelian Higgs model in four dimensions. 
It is described by the action |T^ 

S = f d^Xy/^C (1) 



C = l-D,<f>D'^<f>* + - 7^) - ^(0*0 - v')' + -^R (2) 

Z 4 iOTTCjr 



with 



7^ 



\ 2/32 

where -D^ = V^ — ieA^ is the gauge covariant derivative, is the gauge potential 
of the U(l) gauge symmetry, F^^, is the corresponding field strength and is a 
complex scalar field with vacuum expectation value v. The Born-Infeld coupling 
constant is noted (3^. The geometry is introduced as usual by means of the Ricci 
scalar R. We use the same notations as in 

In the following we study the classical equations of the above field theory in 
the static case and within the cylindrically symmetric ansatz. The metric and 
matter fields are parametrized in terms of four functions of the cylindrical radial 
variable r : 

ds^ = N\r)dt' - dr^ - L\r)dip^ - N^{r)dz^ (4) 
<P = vf{r)e'^^ (5) 
Af.dx^' = ^{n - P{r))d^ (6) 



2 



Here N, L, P, / are the radial functions of r, n is an integer indexing the vorticity 
of the Higgs field around the 2;-axis. 

In the process of establishing the classical equations in terms of the radial 
functions it appears convenient to define dimensionless coupling constants : 

« = y , 7 = 87rGv' , Z^' = ^ , (7) 
and to rescale the radial variable r and the function L{r) according to 



X 



VXv^r , L{x) = VMj^L{r) (8) 



With these redefinitions (and in the units h = c = 1) the diagonal components 
of the energy-momentum tensor read 

= Xv\es + + + u) (9) 

T;*" = \v^{-es + - + + u) (10) 

= Xv'^{es + e^j - iy - e^j + u) (11) 

V = V (12) 

where the different densities are expressed in the dimensionless variable and func- 
tions : 

£. = i(/T , E» = ^ . n^ld-rf (13) 



= . f„ = lg^ . n^^i^^ (14) 

the prime indicates the derivative with respect to x. The corresponding compo- 
nents of the Ricci tensor are given explicitely e.g. in |T3| . 



Then, after an algebra, the classical equations corresponding to the above 
lagrangian and ansatz can be shown to be consistent and lead to the following 
system of four equations for the functions A^, L,P, f : 

(LNN'y f., , P'2 1 1, 
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J^(^)' = afP (17) 

The problem posed by these differential equations is further specified by the 
following set of boundary conditions 

N{0) = 1 , A^'(O) = (19) 

L(0) = , L'(0) = 1 (20) 

P(0) = n , lim P(x) = (21) 

/(O) = , jim/(x) = l (22) 

which are necessary to guarantee the regularity of the configuration on the sym- 
metry axis and the finiteness of the inertial mass defined next. 

The classical solutions can be characterized, namely, by their magnetic field 

eL[r) ar a ovrG L[x) 

and by their inertial mass and Tolman mass by unit length |T^ 

Min = I dXidX2 (24) 

/ iV - V - V - ^^i^^2 (25) 



M 



to 



which in turn can be reduced to one dimensional integrals : 

GM,n = Im^u (26) 

o 

= ^ dxNL{es + e^ + e^ + u). (27) 

GMt, = |m,„ (28) 

= - dxN^L{i^~e^-u). (29) 
8 JO 

In the flat case (7 = 0) the classical energy E of the solutions is obtained directly 
in terms of the inertial mass : E = (27ri;^)Mj„. 
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3 Some special cases. 



In this section we briefly review a few particular solutions of the system above 
and comment on their main properties. 

3.1 Nielsen- Olesen vortices. 

We first consider the equations for the case 7 = 0, /5 = oo. This case corresponds 
to the flat Maxwell- Higgs system; Eqs. (|T5|) , (|l^) are solved by N{x) = 1, L{x) = 



X. The solutions of the remaining equations are embeddings of the celebrated 
two-dimensional Nielsen-Olesen vortices into the three-dimensional space. 
They exist for all values of a and n and possess the following behaviour around 
the origin 

P{x) = n - + 0(a;2+2-) , f(x) = box'' + 0(x"+') (30) 

and asymptotically 

P{x) = a^e--^ , f{x) = 1 - /ooe-"^ . (31) 

The parameters ao, &0; o-oo-, foo are determined numerically. For fixed n, the classi- 
cal energy, say E{a, n), of the solution decreases monotonically when a increases. 



The value a = 2 is very particular because first order equations exist IjT^ which 
imply (|17D, (|18D (they are usually called self dual equations). 

In relation to the self dual case, a = 2, the Nielsen-Olesen solutions enjoy the 
following remarkable property : 

-E(a,n) > E(a,l) for a<2 
n 

^E{a,n) = E{a,l) for a = 2 

-E{a,n) < E{a,l) for a>2 (32) 
n 



as illustrated by Fig. 1 for n = 1,2. This supports the idea |T5| that, in the 



domain a > 2, the n-vortice solutions (with n > 2) are stable against decay into 



n separated solutions with unit vorticity (remembering the exponential decay 
(|3lD of the solutions under consideration) . 



3.2 The Melvin solution. 

Cancelling the Higgs field and the Higgs potential in the initial lagrangian, we 
are left with the Einstein-Maxwell lagrangian. The corresponding equations can 
be obtained from (p!5|)-(|T^) after an appropriate rescaling of x and L : 



X , L \ —L (33) 

V a 

and setting 7 = afterwards (then Eq. ([I8| ) decouples) . A solution of these equa- 
tion was discovered by Melvin |T^. It is namely characterized by the following 
(asymptotic) power law behaviours of the various fields 

N{x 00) = Ax^ , L{x ^ 00) = Bx~^ , P{x 00) = Cx~^ (34) 

where A, B, C are constants (the numerical evaluation gives A ^ 2.25, B ^ 0.44). 
The corresponding masses are Mj„ = Mjp = 2. 



3.3 Gravitating vortices 

The solutions corresponding to the Maxwell case (i.e. = 00) were studied 
at length in [|13|, [l^. We just mention that global solutions (i.e. defined for 
X G [0,00]) exist only on a finite domain in the (quarter of) plane 0,7. More 
precisely they exist for 

0<7<7e(a,n) (35) 

The function 7c(a,n) increases monotonically with a and, around the self-dual 
value a = 2 it roughly behaves accordingly to 7c(a, n) = 2/n+/t(a— 2) , k, ^ 0.15. 
In fact, for fixed a, there exist two branches of solutions for the set of 7 in the 
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domain (|35|). Following |T3[ we refer to the two branches as to the string and the 
Melvin branches. 

For 7 > 7c(a,n), solutions still do exist but they are closed [|l^]; i.e. they 
are only defined for x G [0,Xmaa]- We will not be concerned with these type of 
solutions in the present paper. 



4 Numerical results. 

We solved numerically the system (p!5D-(p!8D for numerous values of the coupling 
constants a,/?, 7, studying in particular the response of the three types of so- 
lutions mentioned above to the Born-Infeld interaction. We now discuss then 
results in the same order as before. 



4.1 The Nielsen-Olesen-Born-Infeld vortices 

We start by considering 7 = 0. The formula ( pT]) still holds for (3 < 00 and, as 
expected, the Nielsen-Olesen solutions are smoothly deformed by the Born-Infeld 
term. Fig. 2 shows the classical energy by means of Mj„(a, (3,n)/n a,s a, function 
of for several values of a and n. It clearly indicates that the inertial mass 
depends only little on (3. The corresponding values in the Maxwell limit 
= 00) are indicated nearby. The main feature demonstrated by this figure is 
that the solutions exist only for a finite range of the parameter /?, i.e. for 

I3c{a,n) < 13 < 00 (36) 

confirming the result of |^. An in-depth analysis about the reasons for this 
critical phenomenon is provided by Fig. 3 where the numerical evaluation of the 
parameters ao, &o for the solutions are plotted as functions of (3'^. Our numerical 
analysis strongly indicates that 

lim an = 00 (37) 
7 



so that, very likely, no solutions exist for (3 < (3c- This critical behaviour is, 
however, not reflected by the value of the classical energy and by the parameter 
bo which remain finite for (3 —>■ j3c- 

Physically, this phenomenon is reflected by the magnetic field of the solution. 



Indeed, from the definition (^31) , it is immediate to see that -B(O) (i.e. the value of 
the magnetic field on the symmetry axis of the solution) is directly related to oq; 
Eq. (|37|) can therefore be rephased in saying that the magnetic field on the z-axis 
becomes infinite when the critical value of (3 is approached. This phenomenon is 
further illustrated by Fig. 4 where the magnetic function B{x) is compared for 
several values of (3"^ in the case a = 1 , n = 1 . 

Refering to the relations (P^j), a natural question is to analyse how the domain 
of stability of the n > 2-vortices (i.e. 2 < a < oo) develops with the Born-Infeld 
interaction. An important element of answer is given by Fig. 1. Clearly the effect 
of the Born-Infeld term is that the stable phase appears for lower values of a. In 

2 

terms of the Higgs particle mass Mfj = \v (remember a = ^) this means that 
the stable phase persists up to a higher value of the Higgs particle mass. 

This result motivates an analysis of the effect of the Born-Infeld term in the 
Georgi-Glashow model for multi-monopole solutions. It is known that in this 
model, self dual equations exist for A = (i.e. for a massless Higgs particle) and 
that multimonopoles solutions are unstable for A > ]TB|. 

It is likely that the Born-Infeld interaction generates a small domain of the 
Higgs mass parameter where the mass of the two-monopole is lower than twice the 
mass of the single monopole, thus leading to a stable phase. Such a phenomenon 
was observed recently in the gravitating Georgi-Glashow model [0. Due to 
gravity, a region of the parameter space exists where the multimonopoles seem 
to be stable against decay into monopoles. 
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4.2 The Melvin-Born-Infeld Solution. 

The Melvin solution gets smoothly deformed by the Born-Infeld term. Decreasing 
we observed that the solution disappears for (3l ~ 6.4; again the phenomenon 
is related to the fact that the value -B(O) tends to infinity for (31. The 

profiles of N{x), L{x) and of the magnetic field for the Melvin solution and its 
deformation for = 7 are compared on Fig. 5 (the dotted line for B has 
B{0) ^ 40). 

4.3 The Born-Infeld gravitating strings. 

The study of the equations (^^-|T^ for generic values of a, /5, 7 is a vast task. Here, 
we have attempted to confirm the features which occured for the flat case. We 
limit ourselves to the case a = 2 but we expect that the results are qualitatively 
the same for the other values of this coupling constant. Our numerical analysis of 
the equations conflrms that, for 7 flxed, the solutions on the string branch as well 
as the ones on the Melvin branch terminate at flnite values of Fig. 6 shows 
the domains of existence of the two types of solutions in the /?^,7 plane. The 
flgure clearly reveals that the value (3^, where the string branch terminates (the 
lower curve) varies only little with 7 and that is is smaller than its counterpart, 
say Pi, for the Melvin branch (upper curve). 
Taking e.g. a = 2, 7 = 1 we flnd 

Pl^OA , 42^5.1 (38) 

respectively for the solution of the string and of the Melvin branches. The inertial 
mass deviates only slightly from the exact value Mj„ = 2 in the Maxwell limit, the 
Tolman mass decreases slowly from ~ 0.62 in the Maxwell limit to ~ 0.4 

for (3' = ~Pl 

For both branches, the value -B(O) of the magnetic fleld on the z-axis becomes 
very large (probably inflnite but our numerical analysis does not allow us to be 
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absolutely sure) when the critical value of is approached. Other quantities, 
e.g. the inertial and Tolman masses remain finite. 

We finally mention that the string branch (which exists for < 7 < 2 in the 
case a = 2) extends slightly into the region 7 > 2 when the Born-Infeld term is 
present. 



5 Final remarks 



To finish, we would like to present a generalization of the result in to the 
gravitating case. The self-dual potential and the corresponding equations (i.e. 
setting A = e^/2 in Eq.(|^) and a = 2 in Eqs.(|T5|)- ([l6|) ) can be replaced by a 
/3^-depending one for which first order equations exist which imply the classical 
(second order) axially-symmetric equations. In the reduced variables, this can be 
achieved by the following substitution 



Vif) = \if - ir —> VsM) =P'\^l-^l- ^(/2 - lyj (39) 
With this potential the self dual equations read 

N{x) = l , L' = l + liPil-f) + e) , Y = J^ ' ^' = ^ 

and imply the full axially symmetric equations of the theory. Further exact 
solutions and algebraic properties can be obtained with this modified theory. 
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Figures Captions 



Figure 1 

The ratio E/n is plotted as a function of a for the Maxwell and Born-Infeld 
vortices and for n = 1, 2. 

Figure 2 

The ratio E/n for the Born-Infeld vortices is plotted as a function of for 
different values of a and n. On each branch, E^ gives the energy of the 
Maxwell solution. 

Figure 3 

The parameters Oq, bo entering in the solutions via Eq. ( pO] ) are reported as 
functions of for different values of a and of n. 

Figure 4 

The profiles of the magnetic field B{x) is presented for several values of (3'^ 
approaching the critical value ~ 0.318. 

Figure 5 

The profiles of the functions L, N, B for Melvin solution (solid lines) and 
Melvin-Born-Infeld solution (dotted lines). On the dotted line -B(O) ~ 42. 

Figure 6 

The domain of the gravitating string and Melvin branches of solutions for 
a = 2. 
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Fig. 2 
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